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Abstract. We investigate the spatial distribution of temperature induced by a 
dc current in a two-dimensional electron gas (2DEG) subjected to a perpendicular 
magnetic field. We numerically calculate the distributions of the electrostatic potential 
</> and the temperature T in a 2DEG enclosed in a square area surrounded by insulated- 
adiabatic (top and bottom) and isopotential-isothermal (left and right) boundaries 
(with 0i o ft < </> r ight and Xi c f t = T r i g ht), using a pair of nonlinear Poisson equations (for 
</> and T) that fully take into account thermoelectric and thermomagnetic phenomena, 
including the Hall, Nernst, Ettingshausen, and Righi-Leduc effects. We find that, in 
the vicinity of the left-bottom corner, the temperature becomes lower than the fixed 
boundary temperature, contrary to the naive expectation that the temperature is raised 
by the prevalent Joule heating effect. The cooling is attributed to the Ettingshausen 
effect at the bottom adiabatic boundary, which pumps up the heat away from the 
bottom boundary. In order to keep the adiabatic condition, downward temperature 
gradient, hence the cooled area, is developed near the boundary, with the resulting 
thermal diffusion compensating the upward heat current due to the Ettingshausen 
effect. 
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1. Introduction 



Thermoelectric and thermomagnetic phenomena [1| have recently been attracting 
renewed interest not only as a route for potentially highly efficient device application, 
e.g., in refrigeration or generating electricity, but also as an effective tool to explore 
fundamental properties of solid state materials j^-H]. Being sensitive to the energy 
derivative of the electric conductivity (or the density of states) or to the entropy of the 
system, thermoelectric and thermomagnetic properties provide us with information on 
the materials complementary to, and often with higher sensitivity than, the information 
inferred from the electric conductivity (215]. For instance, it has been shown that the 
Seebeck or Nernst coefficient measured in a bismuth single crystal 13|, ll4J or a two- 
dimensional electron gas (2DEG) 15[ exhibits clearer quantum oscillations due to the 
Landau quantization compared to those of the electric conductivity (the Shubnikov-de 
Haas oscillations). Moreover, thermoelectric effects introduce additional twist to the 
measurement of electric conductivity or resistivity, or more generally to the distribution 
of the electrostatic potential and electric current. For example, it is necessary to take 



into consideration the thermovoltages in the precision resistivity measurement 16 



Further complication arises by the application of a magnetic field |17l-l20l]. Nontrivial 
distributions of the potential and current, and hence the temperature, can be generated 
by the thermoelectric and thermomagnetic effects. 



Fujita et al. 2l| recently reported seemingly anomalous behavior of the Nernst 



signal in a quantum Hall system, which suggests possible cooling of the electron 
temperature by the current intended to heat the electron system and introduce the 
temperature gradient. Figure [1] (a) shows a schematic diagram of the experimental 
device fabricated from a conventional GaAs/AlGaAs 2DEG. The top (horizontal) bar is 
used as a heater; Joule heating by the heating current 1^ = 4 nA//zm raises the electron 
temperature T ej hi g h there and introduces temperature gradient toward the Ohmic contact 
pad below. The pad is thermally connected to the mixing chamber of the dilution fridge 
kept at T ei i ow = 40 mK, in which the sample is immersed. Thermoelectric voltages are 
measured in the main (vertical) Hall bar. Arms to measure the transverse thermoelectric 
(Nernst) voltage V yx are shown in the figure. With this current heating technique, one 
can heat up the electron temperature selectively, leaving the lattice temperature intact 
(so long as Jh is kept low enough ). T herefore, one can pick out the diffusion contribution 
in the thermoelectric voltages 22| and can eliminate the phonon-drag contribution, 
which is often the dominant contribution in the thermoelectric powers in a 2DEG j^]. 
(Note, however, that the phonon-drag contribution can be negligibly small at T = 40 
mK anyway, see below). The (negative) gate voltage V g applied to the front gate (shown 
by the gray square) allows us to control the carrier density and hence the resistance of 
the heater section independently from the main Hall bar. Figure []Jc) shows the Nernst 
voltage V yx plotted in the B-V g plane. Thermoelectric voltages vanish when the 2DEG 



is in the quantum Hall states [5|, |23|, |24| and Joule heating does not work when the 
heater section is in the dissipationless state. Therefore, nonvanishing signal appears 
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only when both the main Hall bar and the heater section are in between two adjacent 
quantum Hall states, namely only when the Fermi energies (Ep) of both sections cross 
the (disorder-broadened) Landau levels having finite densitiy-of-states. 

Below 5 = 1.8 T, the Nernst voltage V yx behaves as expected, showing oscillations 
as a function of B (taking negative then positive values when Ep of the main Hall 

but does 
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25 



bar crosses a Landau level with the increase of the magnetic field) 
not depend much on V g insofar as the heater section is in the dissipative states. Note 
that V g alters only the heater section and should have no effect on the main Hall bar 
(the section where the thermoelectric voltages are measured). Anomalous behaviors 
are observed above 1.8 T: V yx alternates the sign when V g is swept at a fixed B (see 
figure QJd)). Negative V yx appears for smaller (more negative) V g at the higher magnetic 
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Figure 1. (a) Schematic drawing of the experimental device 2l|. A magnetic field B is 
applied perpendicular to the 2DEG plane. The Nernst voltage V yx is measured between 
the probes depicted as the black squares. The gray square indicates the front gate. 
The current Zj, — 4 nA//im heats the region underneath the front gate to T 6i hi g h and 
introduces temperature gradient toward the bottom pad held at T e ,low = 40 mK. (b) 
Geometry used in the calculation to approximate the heater section (section beneath 
the front gate). Top and bottom (blue) boundaries are assumed to be insulated 
and adiabatic, and left and right (red) edges have fixed potential and temperature, 
0icft = nV, bright = 80 nV, Ti c f t = T r ight = 40 mK. We introduce the polar 
coordinate (r, (p) with the origin located at the left-bottom corner, (c) Color plot 
of experimentally obtained Nernst voltage V yx in the plane of the magnetic field B 
and the applied gate voltage V g , with blue and red colors representing the positive and 
negative values, respectively, (d) Cross sections indicated by solid and dashed vertical 
lines in (c). Anomalous behavior discussed in the text is highlighted by dotted ellipses 
in (c) and (d). 
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field side where V yx is expected to be positive (the areas indicated by dotted ellipses 
in figures [Tfc) and (d)). The inversion of the sign of V yx can naively be interpreted as 
resulting from the inversion of the sign of the temperature gradient, implying that the 
"heater" section can actually be cooled by 1^, depending on the value of V g . 

This speculation led us to investigate, in the present paper, the current- induced 
temperature distribution of a 2DEG placed in a magnetic field, in pursuit of the 
possibility of the current-induced cooling. We numerically examine the spatial 
distribution of the temperature T in a rectangular segment of 2DEG that simulates the 
heater section of the experimental device, fully taking into account the thermoelectric 
and thermomagnetic effects. In order to focus on the current-induced cooling, however, 
we consider simplest possible settings and do not make an attempt to reproduce in 
detail the device configuration (figure Ufa)) of the experiment delineated above; we 
leave out the main Hall bar altogether and approximate the heater section as a square 
surrounded by insulated-adiabatic (top and bottom) and isopotential-isothermal (left 
and right) boundaries in the calculation, as depicted in figure[]Jb). For the method of the 
calculation, we basically follow the prescription developed by Okumura and coauthors 
261 ] for 3D semiconductors at the room temperature and applied their treatment to 
a 2DEG at a low temperature. In this treatment, the contributions of phonons are 
neglected, which is justified at the extremely low temperature (40 mK) considered in 
the present paper. In fact, the electron-phonon scattering time r e _ p h in a GaAs/AlGaAs 
heterostructure in which the 2DEG resides is estimated to be ~ 0.01 — 10 ms at T = 40 
mK 27- 29j]. Despite the wide variation among the literatures, r e _ p h is still orders of 
magnitude larger than the elastic scattering time due to the impurity scattering (or the 
boundary scattering in the case of a small sample in a ballistic regime). The distributions 
of the electrostatic potential (with the electrochemical potential given by — e<ft] see 
the discussion below) and the temperature T are obtained by solving the nonlinear 
Poisson equations, V 2 = F(T, VT, V0) and V 2 T = G(T, VT, V4>), with the functions 
G and F derived from the transport equations [l|, 0, 30-32], as will be detailed in 
section [2j We find that a magnetic field distorts equipotential lines and generates an 
uneven temperature distribution with high- and low-temperature areas emerging at 
the opposite corners of the square (see figure [3]). The low-temperature area is found 
to become colder than the isothermal boundaries, as briefly reported in our previous 
publications 33j,|34|. Two opposite corners with the temperatures lower and higher than 
the bath, respectively, resembling those presented here was previously also reported in 



numerical calculations by Ise et al. [19| when the current is small enough, although the 
mechanism responsible for the cooling was not explicitly specified in that paper. The 
emergence of the cooled part possibly give a qualitative account of the experimentally 
observed sign reversal in the Nernst signal 21[ mentioned above. In order to clarify the 
relation of our numerical results to the experiment described above performed in the 
quantum Hall regime, however, much improvement has to be made in our theoretical 
treatment, as will be discussed in section HI 

The main purpose of the present paper is to clarify the mechanism of the cooling 
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phenomenon found in our numerical calculation. We evaluate the terms in the right hand 
side of the nonlinear Poisson equations F(T, VT, V0), G(T, VT, V0), and identify the 
dominant terms that induce the cooling. We find that the cooling is mainly attributable 
to the adiabatic condition at the bottom edge. It causes a temperature gradient and 
hence the thermal diffusion to balance out the heat current away from the edge driven 
by the Ettingshausen effect, and consequently generates the cooled area adjacent to 
the edge. We will also show with simple analytical arguments that the cooled area, in 
principle, emerges by the application of an arbitrarily small magnetic field. 

The paper is organized as follows. In section [2J we describe the method of 
calculating the spatial distributions of <fi and T by solving the nonlinear Poisson 
equations followed by the results of the calculations. In section [31 we pinpoint the 
mechanism of the partial cooling through the identification of the dominant terms in 
the equations. In section HI we discuss our results in connection with the experiment 
that motivated our study. Section [5] is devoted to conclusions. 



2. Numerical Calculation 



The transport equations describing the electric current density J and the thermal 



current density Jq for isotropic systems [H, 0, 30-32 



arc: 



— V0 = pJ + R(B x J) + aVT + N(B x VT), (1) 
J Q = aTJ + NT(B x J) - kVT + kM(B x VT), (2) 

where B denotes the magnetic field, p the electric resistivity, R the Hall coefficient, a 
the Seebeck coefficient, N the Nernst coefficient, k the thermal conductivity, and M the 
Righi-Leduc coefficient. The transport coefficients, p, R, a, N, k, and M are all defined 
in the isothermal conditions. The terms on the right-hand sides of equations (CEJ) and (J2J) 
respectively represent transport phenomena as follows: Ohm's law, the Hall effect, the 
Seebeck effect, and the Nernst effect in (JTJ; the Peltier effect, the Ettingshausen effect, 
Fourier's law of the thermal conductivity, and the Righi-Leduc effect in (J2]). 
We define the energy-flux density Ju as 

Ju = J Q + 4>J- (3) 

Here we selected the Fermi level (chemical potential at T = 0) as the origin of the energy. 
Noting that the temperature dependence of the chemical potential is negligibly small 
at the low temperatures considered in the present paper, the electrochemical potential 
is given by — e<p, hence the definition ([3]). From equations ([I])-® and the equations 
of continuity in the steady state V • J = V • Ju = 0, we obtain the nonlinear Poisson 
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equations (see Appendix A for the derivation) 
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^G(T,VT,V0). (6) 

The electric current J = (J x , J y ) to be replaced in equations fl3| and (jSJ) is obtained by 
inverting (TXT) : 



J y p 2 + R 2 B 2 




-d x (j) — ad x T + NBdyT 
-d v d> - ad v T - NBd r T 



(7) 



As illustrated in figure (U^b), the left and the right edges are isopotential and isothermal, 
i.e., and T are fixed (the Dirichlet conditions). On the top and the bottom edges, we 
set insulated and adiabatic conditions; the normal components of J and Jq vanish at 
these boundaries. Let us denote the quantities at the top and bottom boundaries with 
tilde, 4>, J, T, and Jq. Equations ((H) and (121) at these boundaries are then reduced to 

-dj = pJ x + ad x f - NBd y f, (8a) 
-d y 4> = RBJ X + ad y f + NBd x f, (8b) 
=NTBJ X - ndyf + nMBd x f. (8c) 
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Figure 2. Discretized two-dimensional sample, where the sample size is L x = L y = 
10 fim and the mesh size is d = 0.1 /im. Thus, the number of the grid points 
is (N x + 1) x (JVj, + 1) = 101 x 101. The label (i, j) denotes the (i, j)th grid 
point along the x- and y-axes, respectively. The downward (blue) and upward (red) 
triangles on i = —1 and on i = N x + 1 have potentials and temperatures fixed to 
OAleft) Tieft) and (bright, Tright)) respectively, set as <^ left = nV and <fi light = 80 nV, 
T\ e ft — T r jght = 40 mK. The (green) squares on the top and bottom boundaries are 
insulated and adiabatic. The magnetic field B is applied perpendicular to the sample. 



We thus obtain the following derivatives (the Neumann conditions), 

a NT ~ 

d y (f) = -(R+ —)BJ X - (oM + N)Bd x T, (9a) 
d y f =^-^-J x + MBd x f, (9b) 

Kj 

where the electric current J x on the boundary is given by 



K 



pn - N 2 B 2 T 



d x (f)+(a-MNB 2 )d x T . (9c 



We numerically solved the set of equations ([!])-© self-consistently on a discretized 
sample illustrated in figure [2] with the Dirichlet and Neumann conditions mentioned 



above. We employed the successive over-relaxation (SOR) method [35|, |36[ in the finite- 
difference calculations. We used the following strategy to achieve efficient convergence in 
the calculation. We started our calculation without a magnetic field, giving the initial 
distributions of cf) and T as linear functions between the left and right isopotential- 
isothermal boundaries. Once we reached the convergence, we then increased B from 
to 10~ 3 T and sought the convergence. We successively increased B step by step with 
an increment of 10~ 3 T, using the result for the previous B as the initial distribution 
for the next B. Nevertheless, we have so far reached only up to 0.2 T, the convergence 
becoming increasingly slower with increasing magnetic field. 
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Table 1. Parameters used for the calculation. 


Parameter 


Value 


Parameter 


Value 


Resistivity p 


20.0 fl a 


dp/dT 


1 nftK- lC 


Hall coefficient R 


-1600.0 flT- ia 


dR/dT 


1 nOT^K- 10 


Seebeck coefficient a 


-0.175 fiVK- lh 


da/dT 


-4.38 pVK- 2h 


Nernst coefficient N 


-7.35 /iVK- 1 T" lb 


dN/dT 


-184 ^VK- 2 T- lb 


Thermal conductivity k 


42.7 pWK~ lb 


dn/dT 


1 nWK- 2b 


Righi-Leduc coefficient M 


-70.0 m 2 V- 1 s- lb 


dM/dT 


b 



a Experimental parameters taken from [21| . 

b Values calculated with semiclassical theories at T = 40 mK using experimentally 
obtained p and R, assuming that the scattering time r depends on the energy e as 



c Tentative values for the numerical calculation, assumed to be small enough not to 
affect the result of the calculation. 



We performed the calculations with the parameters listed in table [TJ p and R are 



taken from the experimental data [21| and other parameters, a, N, k, M, and their 
temperature derivatives are calculated by substituting the values of p and R into the 
semiclassical formulas [5] (for a and N) and making use of the Wiedemann- Franz law 



37| (for k and M) valid for the diffusion contribution. As an initial step, we used the 



values of the parameters at B — throughout the calculation in the present study, 
neglecting their magnetic-field dependence. Although this appears to be very crude 
approximation, the parameters at B = represent roughly the right order of magnitude 
for their values under finite B when the Fermi energy Ep crosses the Landau levels (by 
contrast, some of the parameters vanish in the quantum Hall states), and we believe 
they suffice to discuss qualitatively what happens when the heater section is in the 
dissipative regime. 

Figure [3] shows the distributions of the electrostatic potential 0, the electric current 
J, the temperature T, and the heat current Jq in the case 0i e f t = nV, r i g ht = 80 nV, 
and TJeft = T r ; g ht = 40 mK. We briefly reported the numerical results shown here 



in previous articles 33|, |34|. At B — (figure 0(a)) the electric current J flows 
homogeneously perpendicular to the 0-contours. The temperature, raised by the Joule 
heating, has a symmetric distribution decreasing towards the left and right edges held at 
the fixed temperature. The asymmetry of the heat current Jq, with the net flow going 
to the right, arises because of the left-going potential energy flow due to J (the second 
term in ([3])). Once the magnetic field is switched on (figures Elb)-(d)), the 0-contour 
is distorted. The current J flows nearly parallel to the 0-contours, or more precisely, 
approximately at the Hall angle 

(RB\ 
I 5 (10) 
P J 

deflected from the gradient V0. The Hall angle # H equals 76°, 83°, and 86° at 
B = 0.05 T, 0.1 T, and 0.2 T, respectively. The current J is highly concentrated 
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at the right-top and left-bottom corners (see also figure H]). The distributions of 
<p and J are qualitatively the same as well-known distributions calculated without 
taking the thermoelectric and thermomagnetic effects into consideration 38M41| . More 



quantitative comparison will be made in section 13.11 The temperature distribution 
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Figure 3. Distributions of <j>, J, T and Jq in plane obtained from the 

nonlinear Poisson equations (|3J and ([6]) with the boundary conditions 4>\ e ft = 0.0 nV, 
fright = 80 nV and 7] c ft = T r i g ht = 40 mK. The </>-contours are labeled by the values 
in the unit of nV, and the T-contours by the difference from T] e ft = T^ght in the unit 
of /iK. (a) B = T. The potential gradient generates an uniform current distribution 
and the resulting Joule heat raises the temperature in the mid part, (b)-(d) B > T. 
The magnetic field causes the distortion of the equipotential lines and the distribution 
of the electric current J. It also affects the distribution of the temperature T and 
generates the area having the temperature lower than that of the heat baths. The 
cooled area occurs, in principle, with an arbitrarily small magnetic field, although it is 
apparent above 0.03 T in our discretized sample. (See section f3. 31 for discussion). 
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becomes asymmetric, with high and low-temperature areas emerging around the right- 
top and left-bottom corners, respectively (where J has a high concentration). The 
low-temperature part has temperatures lower than the temperature Ti e f t = T rig ht of the 
heat baths. We will show in section 1331 that, in principle, the cooling appears with an 
arbitrarily small magnetic field. 

3. Mechanism of the cooling phenomenon 

In the present section, we investigate the origin of the cooling phenomenon. 



3.1. Simplification of the governing equations 

As an initial step toward the understanding of the cooling mechanism, we deduce an 
approximate version of the governing equations and the boundary conditions much 
simpler than the original ones. To this end, we compare the terms in the relevant 




Figure 4. (a) Distribution of J. (b) Schematic view of J on the boundaries, (c) 
Distribution of J 2 , (d) Contour lines of J 2 in units of (nA) 2 . All data are for B = 0.1 T. 
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equations using the numerical solutions presented in the previous section, and eliminate 
the terms whose contributions are negligibly small compared to the other terms. First, 
in equation (dJ, we find that |aVT|/|V0| < 1(T 5 and \N{B x VT)|/|V0| < 1(T 6 in the 
magnetic field range (0 — 0.2 T) examined in the present paper. Therefore the terms 
aVT and N(B x VT) can safely be neglected. Noting that dp/dT and dR/dT are also 
negligibly small in a 2DEG at T < 0.1 K, we arrive at the Laplace equation, 

V 2 </> = 0, (11) 
for the electrostatic potential, and the expression 
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p 2 + R?B 2 \ -RB 




(12) 



for the current, as simplified equations to take the place of equations (@J and ([IJ), 
respectively. We can thus calculate with high accuracy the distributions of <fi and J 
neglecting the thermoelectric and thermomagnetic effects. This is to be expected since 
we adopted rather large potential difference but no difference in the temperature between 
left and right edges as the boundary condition. Since the Laplace equation (TTTT) can be 
solved analytical for our boundary conditions H (see seetionEl, this approximation 
vastly simplifies the calculation. 

Next we examine equation ([6]). We find that the term including pj 2 (corresponding 
to the Joule heating) is by far the dominant term, exceeding the other terms by factor 
10 6 . Along with N 2 B 2 T/pn ~ O(10~ 16 ), we have 

V 2 T = -^J 2 , (13) 

K 

as a simplified approximate nonlinear Poisson equation for the temperature. 

The boundary conditions (15a]) - (f5cl) at the top and bottom boundaries can also be 
simplified as, 

-d x (j) = pJ x , (14a) 
-d y ^=RBJ x , (14b) 
=NTBJ X - Kdyf + nMBd x T, (14c) 

(Note that f)14cl) is the same as fl8c|) because all the terms are of comparable order of 
magnitude and therefore cannot be neglected.) Or equivalently, we have the Neumann 
condition 

d y 4> = -RBJ x , (15a) 
d y f = —^-J x + MBd x f, (15b) 

J x = —d x (j). (15c) 
9 
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Equation (1131) reveals that the spatial variation of the temperature in the interior of 
the sample is mainly determined by the Joule heating, with the thermoelectric and 
thermomagnetic effects playing only minor roles. This is not the case at the boundaries, 
where the adiabacity is achieved among the Ettingshausen effect, thermal diffusion, and 
the Righi-Leduc effect, as seen in (114cj) . 

To confirm the appropriateness of the simplification, we calculated <f), J, T, Jq 
with the simplified equations (ITT]) -( ITS"]) and ( TT5]) and obtained the distribution virtually 
indistinguishable from figure [3j The relative differences 5 A = \(A — A')/A\, where A 
and A' are the values from the original and the simplified equations, respectively, were 
sufficiently small for 0, T, J, and J Q as 5(f) < 0(1(T 5 ), 5T < 0(1(T 10 ), 5J X < O(10" 6 ), 
5J y < 0(1(T 6 ), SJ Qx < 0(1(T 7 ), and 5J Qy < 0(1(T 7 ). It verifies that the simplified 
equations (TTTj)-(fl3]) and ( [15]) effectively give the same distributions as those from the 
original equations (SJ)— ^TJ) and (J9j). 

3.2. The role of the Ettingshausen and the Righi-Leduc effects 

In this section, we examine the role played by the first and the third terms in ( I14cp (the 
Ettingshausen and Righi-Leduc effects). We substitute into the coefficients N and/or 
M, and see whether cooled area appears or not. We confirm that it is the Ettingshausen 
effect at the boundary that is indispensable for the cooling effect. 

First, we consider the case N = M = in (I14cl) ; that is, d y T = at the top 
and the bottom boundaries. Figure |5(a) shows the distributions of (f), J, T and Jq 
in this case. Here and in what follows, we used a magnetic field B = 0.1 T as a 
typical example. The result shows that the right-top and the left-bottom corners have 
much higher temperatures than elsewhere owing to the Joule heating by the highly 
concentrated electric-current density J illustrated in figure HI 

Next, we assume that N ^ and M = in fl!4cj> . namely, = NTBJ X — nd y T. 
In this case, we obtain a cooled shown in figure [5(b). It is obvious, therefore, 

that the first term in f ll4cf> (the Ettingshausen effect) plays a major role in the cooling. 
The temperatures on the top (bottom) boundary are now higher (lower) compared 
with those in figure [5(a). This change in the temperature map results from an upward 
temperature gradient d y T = (NTB / k)J x on the adiabatic boundaries (with the negative 
value of N). The temperature gradient can be viewed as being generated so that the 
resulting downward thermal diffusion —nd y T cancels the upward heat current brought 
about by the Ettingshausen effect NTBJ X at the adiabatic boundary, as illustrated in 
figure El If we set, in turn, iV = 0, M ^ in flTIcl) . i.e., = -d y f + MBd x f, the 
cooling effect does not appear as seen in figure |5](c). The boundary condition yields a 
downward gradient d y T = MBd x T, with M < and d x T > at the boundaries. 

Finally, with N ^ and M ^ in (I14c[) . we obtain the distributions shown in 
figure [5(d), which is the same as figure [3(c) but re-presented for comparison. Since the 
term MBd x T (the Righi-Leduc effect) reduces d y T, a cooled area becomes smaller than 
that for N ^ 0, M = (figure [5(b)). 
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The distributions of and J remain virtually unaltered throughout figures Efa)-(d), 
despite the change in the boundary condition f)14cl) . This is because they are basically 
decoupled from the thermoelectric and thermomagnetic effects in the present situation, 
as demonstrated in section 13.11 

3.3. Threshold magnetic field for the cooling 

In this section, we show that a cooled area is generated, in principle, by an arbitrarily 
small magnetic field; we can always find a cooled area if we can approach indefinitely 
close to the left-bottom corner. In practice, however, the minimum distance from the 




(a)jV = Af=0 



20 40 60 80 100 

(b)JV^0,M=0 



(c),V 0, Af^O 



20 40 60 80 100 



(d)JV#0, M^0 



Figure 5. Distributions of 0, J, T and Jq in plane at B = 0.1 T for (a) 

N = M = 0, (b) N ^ 0, M = 0, (c) N = 0, M ± and (d) N ^ 0, M + 0. (The 
data in (d) is the reproduction of figure HJc).) The ^-contours are labeled in the unit of 
nV. The labels for the T-contours indicate the differences from Ti c f t = T r i g ht = 40 mK 
in the unit of fiK. 
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^ «^ 1 i A Ml Hi 




kMB3 x T d x T 



NTBJ, 



-jcd v T 
NTBJ X y 




kMB3 x T d x T 
_tii_ 



Figure 6. Schematic diagram showing the thermal current flow around the left- 
bottom and right-top corners. On the top and bottom boundaries, the total heat 
current, which is the sum of the heat current by the thermal diffusion {—nd y T, red), 
the Ettingshausen (NTBJ X , blue), and the Righi-Leduc (nMBd x T, green) effects, 
vanishes. The heat currents by the Ettingshausen and Righi-Leduc effects go upward 
and downward, respectively, since J x < 0, d x T > 0, N < 0, and M < 0. The right-top 
and left-bottom corners have temperatures higher and lower than the right and left 
boundaries, respectively, since d y T > 0. 



corner is limited by a certain physical length scale (obviously, the length, e.g., much 
smaller than the inter-atomic distance of the host crystal does not make sense), which 
sets a threshold for the magnetic field to generate the cooled area. In our discretized 
system used for the numerical calculation, the minimum distance is the separation 
between the grid points (not the physical length scale but rather an artificial distance), 
which determines the lowest magnetic field for the cooled areas to be observed on our 
grid points. 

As described in section [37T| equations (0|-([7]) and (J9]) are well-approximated by the 
simplified equations (fTTI) (jT3|) and ( fl~5l) for the present system. Equation ( fill under 



the boundary conditions f)15ap and f ll 5c[) can be solved analytically [39|. Although we 
have limited ourselves to a square sample thus far, analytical solutions are given more 
generally for rectangular samples (L x ^ L y in figure [2]). The electric field E = — V</> is 
given by 

E x = — E^e 1 cost?, (16a) 
E y = £ e 7 sin tf, (16b) 

with 

= -Aft V 1 Sinh ^ 2n - 1 - >?r ^ CQS ^ 2U - jj^i (17) 
H ^(2^-l)vr cosh[(2n-l)f] ' 1 ' 

ti = Aft 1 cosh ^ 2n ~ jjl^l sin ^ 2n ~ 

^(2n-l)7r cosh[(2n-l)f] ' 1 ' 

where we introduced normalized coordinates £ = x/L x , r\ = (y—L y /2)/L x and the aspect 
ratio A = L y /L x . The constant E Q is determined by the potential difference between 



15 



the side boundaries A0 = r i g ht — 0Mt, the Hall angle 6*h in ([TO]) (or the magnetic field), 
and the aspect ratio A as 

A(j) 



E ° I(d n ,a)L x 



I(0 H ,A) 



cos 





J ( sech— 

7T 2 



4gH ^sin[(2n-l)<] se( . 



7T 



7T 



(2n-l)-A [>J£ 



(19) 
(20) 
(21) 



with Jo (a;) the Bessel function of order zero. 

From equations f lT2|) and f|T6|) . we have the electric current density, 

cos#h 



Jx — —Er 



-e 7 cos( , i9 — 6>h), 



P 



P 



, cos 6> H 

p 



(22a) 
(22b) 
(22c) 



We obtain the total current Jtot by integrating J x along an arbitrary axis in y-direction, 

r \/2 



r.X/2 

J tot = -L x / J x (f = const.; rj) dr) 

J-X/2 



L X E - 



-A/2 

cos(6> H ; 



P 



(23) 



„1 / oo 

X(0 H ,A) = A / cos ^h^(-i)"" 1 

J V n=l 

cosh [(2n-l)^V] 



(2n - 1) 



-sech 



7T 



(2n-l)^ 



(24) 



For A = 1 (square sample), it can be shown, to an extremely good approximation, that 
K(6n, 1) ~ J(#h, 1), which monotonically decreases with increasing #h from 1 at 6*h = 
to 0.847 at #h = tt/2. Therefore we have 

A0cos(6» H ) 



Jtot (A 



(25) 



We introduce here the polar coordinate (r, <p), where r = 1 x 2 + y 2 and <p = 
tan _1 (?//x), with the origin located at the left-bottom corner of the system as shown 
in figure GJb), noting that J is nearly isotropic in the vicinity of the corner as can be 
seen in figure H](c) and (d). Using the polar coordinate notation for the temperature 
gradient, r r = d r T and = r^d^T, the Poisson equation (|T3|) is written as, 
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Since we find that r r <C in the vicinity of the corner r = in our numerical result, we 
can neglect the first term in (1261) . Thus, we can express the gradient t v at the isothermal 
boundary (if = ir/2) with a small r by integrating ( 1261) as, 

r,(l,r) = r^ =Q -^ r^rJ^. (27) 



o 



Since J is nearly isotropic, we can replace J in the integral by J x at the bottom boundary, 
Jl^o = —J x \ y=0 = J(r) (note that J x < and J y = at the bottom boundary). With 
this approximation, we obtain 

^f.0 =T *u-^ rJ(r) ' (28) 

Noting that M is not essential for the cooling effect (see section l3T2"j) . we assume M = 
in f)14cp for simplicity to obtain 

= dy T \y=o = ~^r J ^- ( 29 ) 

Using in (|25j|, we have 

/ 7T \ NTB T . . p 7T T . . 9 . . 

^ ( 2 ' V = ~ JW ~ k 2 r J(r) ■ (30) 

As illustrated in figure [TJ temperatures lower than T\ e { t = T(ip = vr/2) emerge if 
T (p (n/2,r) > 0, namely if 

B > * 9 rJ(r) - 7r ^ cosg H rc ^) / 31 n 
2(-iV)T rJ(rj " 2 (-iV)T re ' 



^ T 




Figure 7. Schematic illustration of the ^-dependences of (a) t v = T~ x dipT and (b) 
T around the left-bottom corner. The points at <p = (solid circle) and <p — tt/2 
(open circle) correspond to the bottom and the left edges, repectively. T v (<p — 0) > 
and T(ip = tt/2) = 71 e ft have fixed values determined by the boundary conditions. 
t v decreases with increasing ip, following (|26p. If Tu>((f — ft I '2) > 0, T increases with 
increasing ip near the left edge, ensuring the presence of the area colder than Tj B ft. 
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with 7(r) = 7(£ = r/L x , rj = 0), and we used iV < in the derivation. As will be shown 
below, the right-hand side tends to with r — > 0, although the current J(r) diverges 
with r — > (see figure H]). On the bottom boundary, we have from ( jTTl) . 



7(r) 



4# H (— In 

Z7T 



cot 



7rr 



2L q 



n=l 



COS 



(2n - 1)tt- 



(2n- 1)tt 



1 — tanh 



(2n 



Att 



(32) 



The second term in the large round brackets is less than for r — > 0, and thus from 
a, 



t / \ ^ 77i COS # H / 2^ H , 

J(r) < £/ — : — ex P s In 



E, 



P 

COS 6*h 

p 



Ti 



cot 



7rr 



2L, 



7TT 

2L, 



+ 0(r) 



(33) 



Since < 28u/^ < 1, we have rJ(r) — >■ for r — > 0. 

With a fixed B, an area within the distance r from the left-bottom corner becomes 
colder than the isothermal boundary (ip = vr/2) for r satisfying ( 13T|) . Alternatively, for 
fixed r (~ 0), ( )3T|) gives a threshold magnetic field for the position r to become colder 
than the isothermal boundary, (which can be made, in principle, arbitrarily small by 
letting r — > 0). Practical thresholds for the numerical calculation is given by the mesh 
spacing d. By setting r = d = 0.1 fim, we have a threshold magnetic field 0.03 T, which 
is consistent with our numerical results for B = 0.03 T (not shown). 

4. Discussion 



We have shown that the electron temperature can be cooled down by a dc electric 
current with an arbitrarily small magnetic field. The emergence of the cooled area 
is qualitatively consistent with the cooling inferred by the experimentally observed 
sign reversal of the Nernst signal 2l| described in the Introduction. Since the present 
calculation is performed on much simplified setups compared to the actual experimental 
arrangements, however, it is rather difficult to unambiguously relate the present results 
to the experimental findings at this stage. In this section, we point out noticeable 
discrepancies between the setups in our calculation and the experiment, and discuss 
modifications necessary to be made in the theoretical treatment to fill the gap and to 
facilitate more precise and quantitative comparison in the future study. 

First, the calculated temperature decrement is extremely small, order of 10 -2 /iK. 
It is rather unlikely that the effect caused by such a small temperature change can 
be experimentally detected. It will be necessary to alter our boundary conditions to 
reflect the experimental conditions more precisely. Above all, we set the temperature 
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of boundaries immediately to the left and right of the "heater section" fixed, while in 
the experiment the temperature is fixed at the left and right contact pads (hatched 
rectangles in figure Ufa)) separated ~300 //m away from the heater section (the region 
below the front gate depicted by the gray rectangle in figure Ufa)). It seems plausible 
to expect from the mechanism described in section [3] that the separation allows larger 
variation of the temperature at the left-bottom corner of the heater section. 

Our boundary condition also neglects the possible effects brought about by the 
heating current Jh passing through the regions with different carrier densities (regions 
with and without the front gate, see figure EUa)) hence having differing transport 



coefficients. In fact, it has been pointed out 42|, |43] that the heat can be either absorbed 



or emitted by 1^ at the boundary between the regions having different filling factors, 
owing to the difference in the Peltier constant in the two regions, resulting in cooled or 
heated areas. This can be an alternative candidate for the origin of the cooling effect 
distinct from the mechanism discussed in the present paper. 

The fixed potential difference between the right and left boundaries A0 is also at 
variance with the experiment, in which the total current J tot is fixed. According to 
equations ( {TO]) and ( |25|) . J tot decreases with increasing B if A0 is kept constant. We can 
envisage larger temperature change if we increase A0 with B to keep J tot constant as 
in the experiment. (We found difficulty, however, in the convergence of the numerical 
calculation if we set a larger value for A0.) 

So far, our calculation is limited to rather small magnetic field (B < 0.2 T) for 
technical reasons, especially for the bad convergence of the calculation. Noting that both 
the temperature decrement and the area of the cooled region increase with increasing 
magnetic field (see figure 3), we can expect much more pronounced cooling effect if we 
can extend our calculation to larger magnetic field. 

Second, the sign reversal (corresponding to the appearance of the cooled area) was 
observed only when the magnetic field was large enough (B > 1.8 T) and only for the 
value of V g at which electrons occupy less than half of the topmost Landau level (see 
figure 1 (c)). This can be related to the magnitude and the sign of the Ettingshausen 
coefficient NT, or equivalently, to those of the Nernst coefficient N. (Note that the 
Ettingshausen and Nernst coefficients are related by the Kelvin-Onsager relation.) As 
mentioned earlier, we used, in the present calculation, the value of iV(< 0) at B=0 for 
simplicity, neglecting the .B- dependence. In reality, the magnitude \N\ decreases with 
increasing B, until Landau levels are clearly resolved. In a quantizing magnetic field, 
\N\ again takes a large value when the Fermi energy Ep lies in a (disorder broadened) 
Landau level, with the sign of N alternating depending on whether Ep is below or 
above the center of the Landau level, namely, whether the energy derivative of the 
density of states is positive (electron-like) or negative (hole-like) [5J||. It can readily 
be seen from the discussion in section [3] that the cold and hot areas appearing in the 
left-bottom and right-top corners, respectively, interchange their roles when the sign of 

\ Note that N in the present paper corresponds to S xy /B in Q. 
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N is inverted. The conditions for B and V g mentioned above thus can be interpreted as 
the condition that N possesses large enough magnitude (for the Ettingshausen effect to 
have sufficient strength) and the appropriate sign (so that the cooled area is generated 
on the side adjacent to the main Hall bar), respectively. 

Third, in a 2DEG in the transition region between two adjacent quantum Hall 
states, the current is carried by both the bulk extended state of the topmost Landau 
level and the edge states from lower Landau levels. Although the edge states can have 
significant impact on the distributions of the electric and heat currents j-JU^. and hence 
on the temperature distribution, our calculation amounts to neglecting the edge states 
altogether. 

Apparently, much improvement has to be made to quantitatively explain the 
experiment motivated our study. We believe, however, that our simplified approach 
has been advantageous to pinpoint the very essence of the cooling mechanism. 



5. Conclusions 



We have investigated the temperature distribution induced by a dc current in a 2DEG 
subjected to a perpendicular magnetic field, surrounded by isopotential-isothermal 
(left and right) and insulated-adiabatic (top and bottom) boundaries (figure Hfb), 
figure |2]). By numerically solving the nonlinear Poisson equations @ and ([6]), we have 
demonstrated that an area having the temperature lower than the isothermal boundaries 
(kept at 40 mK) appears in the vicinity of one of the corners where the electric current 
density is highly concentrated. The cooling is ascribed to the Ettinghausen effect, 
which pumps the heat away from the adiabatic boundary. The adiabatic condition ( 1Hc|) 
at the boundary requires the temperature gradient to be generated, with the resulting 
thermal diffusion canceling out the thermal current due to the Ettingshausen effect. We 
have shown that, owing to the temperature gradient, the cooled area emerges with an 
arbitrarily small magnetic field, although the area shrinks within closer proximity of the 
corner with decreasing magnetic field. 

Although the present study is motivated by the recent experiment 2l|, the relation 
of the present mechanism to the experiment remains rather unclear. Nevertheless, the 
confirmation of the presence of the rather counterintuitive current-induced cooling effect, 
as well as the identification of the mechanism responsible for the effect, appears to 
be important in its own right, and underlines the complication brought about by the 
thermoelectric and thermomagnetic effects. 
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Appendix A. Derivation of the nonlinear Poisson equations ((4j) and (|6]) 

In this Appendix, we present, for completeness, the derivation of the nonlinear Poisson 
equations (@| and (jSJ) from the transport equations ((Tj) and @ in the steady-state. We 



basically follow the prescription presented in [26[ and [45[. See also [1| 



In the steady state (d/dt = 0), we have, from the continuity equations, 

V-J = 0, (A.l) 

and 

V • J v = V • ( J Q + <j>J) = 0. (A.2) 
Taking the divergence of (pQ), we have, with the aid of (lA.lj) : 



- V 2 =V • [pJ + RB x J + aVT + NB x VT] 
=Vp ■ J + VR ■ (B x J) + R[V ■ (B x J)] 
+ Va • VT + aV 2 T 

+ VN-(Bx VT) + N[V • (B x VT)] . (A.3) 

Noting that, under the temperature gradient VT, an arbitrary transport coefficient 
X depends on (x, y) through its temperature dependence, the gradient of X can be 
expressed as 

dX 

VX = — VT. (A.4) 



We thus can rewrite flA.3j) as follows 



+ — [VT-(Bx J)}+R[V-(Bx J)] 

+ ^(VT) 2 + «V 2 T 
dN 

+ — [VT-{Bx VT)] + N[V-(Bx VT)]. (A.5) 
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Using formulas for vector operation, (1A.5j) becomes 



M {J . VT) + ™ lVT . (BxJ)] 

+ R[J ■ (V x B) - B ■ (V x J)] 
+ ^(VT) 2 + aV 2 T 
+ iVVT- (V x B). 



(A.6) 



Using further the Maxwell equation V x B = pj + ep(dE/dt) = in the steady state, 
omitting the negligibly small term pj, (IA.6[) becomes 

-V 2 0=^(J-VT) 
dR, 



+ -[VT.(BxJ)]-i?[B.(VxJ)]. 
+ ^(VT) 2 + aV 2 T. 



(A.7) 



Here we calculate V x J from under the conditions (1A.1I) and (1A.2I) . The 
rotation of ([!]) gives 



= -V x V0 
= V x [p J + RB x J + aVT + NB x VT] 

= ^(VTxJ)+p(VxJ) 

+ (J ■ V)(i?B) - (i?B • V) J - J[V ■ (KB)] 
+ RB(V ■ J) 

+ ^(VxVT) + q(VxVT) 

+ (VT • V)(JVB) - (iVB • V)VT 
- VT[(V • (NB)} + NB(V ■ VT) 

= ^(VTxJ)+p(\7xJ) 



dR , . m . „ 

+ -(J.VT) B + 



ivv 2 r + ^(vr) 2 

(JjJ. 



B. 



(A.8) 



We thus obtain 



V x J = - 

+ 



1 ( dp 



p { dT 
dR 



(VT x J) 



dN 



— (J.VT) + NV 2 T + —(VT) 2 



B 



(A.9) 
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Substituting ( ]A.9j) into fj A. 7[) . we obtain the following equation: 

+ ^[VT.(BxJ)] + ^[B.(VTxJ)] 

RB 2 \dR /T „ m , , - m di\T 
+ — ^(J-VT) + NV 2 T + —(VT) 2 

+ %{VT) 2 + aV 2 T 



dT 

da RB 2 dN 
dT + p dT 



(VT)' 



_i^_f) [V r.(BxJ)] 



+ 



p dT dT 
dp iZS 2 di? 
dT p dT 



(J ■ VT) 



+ [a + 



RNB' 1 



V 2 T. 



(A.10) 



By substituting ([6]) (to be derived below) in ( 1A.10I) . we arrive at equation (J4]). 

Next, we derive the Poisson equation for T, OH]). We obtain the following equation 
from the divergence of J\j using equations (lA.lj) and f ]A.2j) : 

V • (0 J) + V • (aTJ) + V • [NT(B x J)] 

- V • (kVT) + V • [kM(B x VT)] 
= V0 • J + Va • (T J) + aVT • J 

+ (TVN + NVT) • (B x J) + NT[V ■ (B x J)] 

- Vk • VT - kV 2 T + (MVfc + kVM) ■ (B x VT) 

+ kM[V'(BxVT)]=0. (A.ll) 



Using (|A.4h . we have 



v^.j + (r^ + a )(vr. 



dN 



+ [ T— + N)[VT-(Bx J)} + NT[V ■ (B x J)] 



kV 2 T 



+ 



M^V^VT^xVT)] 



+ kM[V ■ (B x VT)] = 0. 



(A.12) 
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Using again formulas for vector operation, we rewrite the above equation as, 



V0-J+(V^ + a) (VT-J) 
+ ( T^ + N) [VT ■ (B x J)} 



dT 

+ iVT[J • (V x B) - B • (V x J)] 
-^(VT) 2 -kV 2 T 

+ kM[VT ■ (V x B)] = 0, (A.13) 
which reduces, with the Maxwell equation V x B = 0, to 

K V 2 T =V0 ■ J + (t^ + a) (VT • J) 



+ + Nj [VT ■ (B x J)] 



<7k 

— NT[B ■ (V x J)] — — (VT) 2 . (A.14) 

We substitute (CQ) into the above equation and obtain 

K V 2 T = -[ P J + R(B x J) + aVT + N(B x VT)] • J 

+ (^H (vt - j) 
+ + n ) [vt ■ (s x j)] 

- NT[B ■ (V x J)] - ^(VT) 2 
= -pJ 2 + T^(VT-J) 



+ (T^ + 2N) [VT ■ (B x J)] 



f7 K 

— NT [B ■ (V x J)] — — (VT) 2 . (A.15) 
By substituting (1A.9I) into ( 1A.15I) . we arrive at fl6]). 
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